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ABSTRACT 

In this paper we present our results of a comprehensive spectroscopic analysis of 
WASP-1, the host star to the exoplanet WASP-lb. We derive T cS = 6110 ± 45 K, 
log g = 4.28 ± 0.15, and [M/H] = 0.23 ± 0.08, and also a high abundance of lithium, 
log n(Li) = 2.91±0.05. These parameters suggests an age for the system of 1-3 Gyr and 
a stellar mass of 1.25-1.35 Mq. This means that WASP-1 has properties very similar to 
those of HD 149026, the host star for the highest density planet yet detected. Moreover, 
their planets orbit at comparable distances and receive comparable irradiating fluxes 
from their host stars. However, despite the similarity of WASP-1 with HD 149026, their 
planets have strongly different densities. This suggests that gas-giant planet density 
is not a simple function of host-star metallicity or of radiation environment at ages of 
~ 2 Gyr. 
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1 INTRODUCTION 

Close-orbiting giant exoplanets that transit their parent 
stars offer unique insights into the interior structure and 
evolution of gas-giant planets. They are the only extrasolar 
planets whose radii and masses can be determined unam- 
biguously from their light curves and radial-velocity varia- 
tions. The discovery of the first such planet, HD 209458b 
(Charbonneau et al. 2000 Henry et al. 20001, established 



its gas-giant nature. Since then, 13 other transiting planets 
have been discovered, revealing a diverse range of densities 
and hence internal compositions. The Saturn-mass planet 
HD 149026b represents one extreme, having a high density 
requiring a substantial rock/ice core. At the other end of the 
density scale one finds HD 209458b, HAT-P-lb flBakos et all 



2007| and WASP-lb ( |Collier Cameron et a!||2007| ). These 



planets have densities so low that even coreless models have 
difficulty reproducing their radii ( Fortney et alpOOrJ Bur- 
rows et al||2006| ). |Guillot et al.| ( |2006| ) pointed out that the 
core masses inferred for then-known transiting planets show 
a weak correlation with the metallicity of the parent star. 



Burrows et al. ( 2006 1 find a similar trend. They also noted 



that WASP-lb, WASP-2b and XO-lb (McCullough et al 



20061 have very similar masses but substantially different 



* Based on observations made with the Nordic Optical Telescope, 
f E-mail: Eric.Stempels@st-andrews.ac.uk 



radii. In the absence of detailed spectroscopic abundance 
analyses for WASP-1 and WASP-2, however, it is difficult 
to determine a likely cause for their inflated radii. 

WASP-lb has, in comparison with other planets, a very 
low density; it is 2.4 times as massive as HD 149026b, but 
also 3.3 times less dense. This poses a significant challenge 



to existing theoretical models. Burrows et al. (20061 find it 
hard to reconcile the planet's large radius R p = 1.44i?j up 



(Charbonneau et al. 2007 Shporer et al. 20071 with evo- 



lutionary tracks even for coreless models with ten-times- 
solar atmospheric abundances and ages greater than 1.5 
Gyr. This leads to a wide range of speculations. Could it 
be possible that the star (and by implication its planet) 
is relatively young, and thus still contracting? Could the 



strong irradiation of the planet inhibit contraction ( Guillot 
fc Showman||200"2 1? Can additional interior opacity due to 
enhanced atmospheric metallicity slow contraction ( |Burrows| 
et al. 20061? Or perhaps (as a counter-argument) can an 
increased interior molecular weight due to enhanced atmo- 
spheric metallicity accelerate contraction? Can tidal heat- 



ing caused by orbital eccentricity ( Bodenheimer, Lin & 



Mardling 2001) or rotational obliquity (Winn & Holman 
2005[ ) provide additional interior pressure support? Given 
so many theoretical candidate mechanisms for inflating a 
planet's radius, comparative studies of the rapidly-growing 
number of stars that host well-characterised transiting plan- 
ets are essential if we are to identify the dominant environ- 
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Figure 1. Radial velocity measurements from the seven obtained 
spectra of WASP-1, compared with the ephemeris determined by 
|Collier Cameron et aL] | |2007| l. 



mental and evolutionary processes that determine a mature 
gas-giant planet's internal structure and outer radius. 

Here we present a detailed analysis of the spectrum 
of the parent star to WASP-lb, WASP-1 (= GSC 02265- 
00107), with the goal of establishing its metallicity and 
evolutionary status. We use as our starting point the pre- 
liminary stellar parameters published by |Collier Cameron] 
|et al.| ( |2007[ ). Although their analysis is roughly similar to 
the one presented here, the spectra they obtained from the 
SOPHIE spectrograph at OHP were compromised by back- 
ground scattered-light contamination and uncertain contin- 
uum normalization, leaving the metallicity unconstrained. 
We also compare the metallicity, effective temperature and 
evolutionary age to those of the very similar host star of 
HD 149028b, and discuss whether the properties of the host 
stars give any useful clues to the strongly-contrasting den- 
sities of their close-orbiting gas-giant planets. 



2 OBSERVATIONS AND DATA REDUCTION 

WASP-1 was observed in early January 2007 at the 2.5m 
Nordic Optical Telescope (NOT), as part of the science ver- 
ification program of the newly commissioned spectrograph 
FIES (Fiber Echelle Spectrograph). We obtained a total of 
seven spectra on six nights, with a spectral resolution of 
R = 47 000 and an exposure time of 20 minutes each, cover- 
ing the wavelength region 4000-7350 A at a signal-to-noise 
of 40-50 in each spectrum. 

The FIES spectrograph is a bench-mounted fiber- 
fed echelle spectrograph, located in its own temperature- 
stabilized building and without any moving parts. With 
such a particular setup one obtains very stable spectra in 
a fixed format for which one of the authors (HCS) devel- 
oped an automated data reduction system. This system uses 
Python and PyRAF to access the echelle reduction rou- 
tines of IRAF. The reduction is fine-tuned for the properties 
of FIES, and performs all necessary reduction steps, such as 
subtraction of biases and scattered light, flat-fielding, or- 
der extraction, normalization (including fringe-correction) 



Table 1. Obtained parameters for WASP-1. See Sect. [5~3] for a 
description of how the uncertainties were derived. 



Parameter 



Value 



Correlated 
uncertainty 





6110 


± 


45 K 


± 50 K 


logs 


4.28 


± 


0.15 


± 0.05 


[M/H] 


0.23 


± 


0.08 


± 0.03 


logn(Li) 


2.91 


± 


0.05 


± 0.05 


[Na/H] 


0.12 


± 


0.08 


± 0.02 


[Si/H] 


0.26 


± 


0.06 


± 0.01 


[Ti/H] 


0.30 


± 


0.16 


± 0.01 


[Fe/H] 


0.26 


± 


0.03 


± 0.02 


[Ni/H] 


0.24 


± 


0.07 


± 0.02 


[Mg/H] 


0.10 


± 


0.05 


± 0.04 


v sin i 


5.79* 


± 


0.35 kms" 1 




«rad 


-13.46 


± 


0.1 kms" 1 





Assuming a fixed value for the macroturbulcncc. 



Otherwise an upper limit (see Sect. 3.3 1 



and wavelength calibration, leaving the observer with fully- 
reduced spectra ready for scientific analysis. 

Using the associated wavelength calibration frames we 
corrected each spectrum for instrumental shifts. We also ap- 
plied the correction for the heliocentric velocity and checked 
and calibrated the spectra against a radial velocity stan- 
dard. For the purpose of determining the stellar parameters 
of WASP-1, we also created one high-quality spectrum with 
excellent signal-to-noise (~ 100) by combining the seven in- 
dividual spectra into one. 



3 PROPERTIES OF WASP-1 
3.1 Radial velocities 

The seven individual spectra we obtained provide an inde- 
pendent check of the sinusoidal radial-velocity variations de- 



tected in WASP-1 by Collier Cameron et al. (20071. In order 



to obtain the best accuracy, we cross-correlated high signal- 
to-noise regions of the spectra on an order-by-order basis. In 
this analysis we excluded orders that contain atmospheric 
features. This allowed us to determine radial velocities with 
an accuracy of ±125 m s -1 . We find that the ephemeris 



determined by Collier Cameron et al. (20071 is compatible 
with our observations (see Fig. |1|. 



3.2 Spectroscopic analysis 

Several detailed spectroscopic investigations of exoplanet 
hosts have been presented in the literature (among others, 
Gonzalez et al.|[200T] |Valenti fc Fischer|[2005] |Santos, Is-| 



raelian, & Mayor 



2004 



Santos et al. 2006). The primary 



difference between these studies is the choice of numerical 
methods and model atmospheres to determine which stellar 
parameters best describe the stellar spectrum. In this paper 
we choose to follow the methodology of one of these studies, 
namely [Valenti fc Fis cher (2005 , hereafter VF05), which was 
also used in the studies of the transiting exoplanet systems 
XO-1 ( |McCullough et al.pOOo] ) and HAT-P-1 (|Bakos et al 



|2007[ ). We make identical assumptions and use of the same 
tools, techniques and grid of model atmosheres, allowing us 
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Figure 2. The above panels show a comparison of the observed spectrum (grey) with the synthetic spectrum based on the recovered 
parameters (solid, thick line). Panel a) shows the temperature-sensitive He? 6563 line, b) the Mgb 5175 triplet, sensitive to logg c) a 
section of the region containing a wealth of metal lines, sensitive to [M/H] and d) the Li I 6708 line. Thes stellar parameters of WASP-1 
were not determined from each panel individually, but from all spectral regions simultaneously (see text). 



to obtain stellar parameters for WASP-1 that can be directly 
compared to the results of VF05. 

We analysed our high-quality spectrum of WASP-1 with 



SME (Spectroscopy Made Easy, see Valenti & Piskunov 
|1996[ ), an IDL-based program that uses synthetic spectra 
and multi-dimensional least-squares minimization to deter- 
mine the best set of stellar parameters (the effective tem- 
perature T e g, the gravity loggr, the metallicity [M/H], the 
projected radial velocity vsini, the systemic radial veloc- 
ity frad, the microturbulence u m i c and the macroturbulence 
D mac ) that describe an observed spectrum. 

Successful calculation of synthetic spectra requires a 
grid of model atmospheres covering the parameter space of 
interest, as well as accurate line lists of atomic transitions. 
In analogy with VF05, we used for our analysis the routine 
for 3-dimensional interpolation on the [Kurucz ( 1993 1 grid 



of LTE model atmospheres. Atomic line data was obtained 



from the VALD database (Piskunov et al. 1995 Kupka et 
|aLl|l999"l ). To improve our ability to model the stellar spec- 
trum in the spectral regions of interest (see below), we ad- 
justed the oscillator strengths and broadening parameters 
for some of the lines in our line lists. This was done by com- 
paring the NSO spectrum of the Sun ( Kurucz et al.|1984 l to 
a synthetic spectrum (using the parameters T e fi = 5770 K, 
logg = 4.44, vs'mi — 1.4 kms -1 , « ra( j = 0.4 kms -1 (grav- 
itational blueshift), v m i c = 0.866 kms -1 and u mac = 3.57 
kms -1 and solar abundances). We also checked our results 
against a high-quality spectrum of Procyon, obtained with 



FIES at the NOT (for stellar parameters, see Fuhrmann et 
aLl[l997l. 



In order to use assumptions identical to those of VF05, 
we decoupled the correlation between microturbulence w m j c 
and metallicity by fixing the value of v mlc to 0.85 kms -1 . 



Similarly we followed their empirical relation for the value 
of the macroturbulence, giving f mac = 4.5 kms -1 for a star 
with T c ff « 6200 K (see also Sect. |3.3| for additional re- 
marks on the choice of macroturbulence). We also released 
the elemental abundances of five elements (Na, Si, Ti, Fe 
and Ni). As an initial guess we used the parameters deter- 
mined by Collier Cameron et al. (20071, as well as [M/H] 
= 0, vsini = 0.5 kms" 1 , and (Fe/H) = -4.50. 

To constrain the full set of stellar parameters we iden- 
tified three important wavelength regions (5160-5190 A, 
6000-6200 A and 6540-6590 A). The first region contains the 
triplet of Mgb lines, primarily sensitive to logg, while the 
second region contains a large number of well-isolated and 
unresolved spectral lines of a range of different elements, 
sensitive to v la _d, vsini, [M/H], and the individual elemen- 
tal abundances in particular. 

The third region covers the broad Ha 6563 line. Al- 
though this region was not included by VF05, we consider 
this region as a strong indicator to constrain effective tem- 
perature. The broadening of the outer wings of this line 
is sensitive to T e g, and largely independent of other stel- 
lar parameters. Our synthetic synthesis calculations use the 
state-of-the-art hydrogen line broadening theo ry of|Barklem,| 
Piskunov, & O'Maral d2000b. iBarklem et al.l d2002b showed 



that this theory can accurately reproduce the Balmer line 
wings of stars with a range of different spectral types, and 
that it produces good agreement with the results of the in- 
frared flux method (IRFM) developed by |Blackwell, Petford, 



& Shallis ( 19801, but with smaller error margin^] They also 



1 While LTE simulations of Balmer line wings are widely used to 
determine stellar temperatures, Barklcm (2007) recently pointed 
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showed that there is a weak dependence on the model at- 
mosphere used. 

During our analysis we encountered difficulties in ob- 
taining a good fit to the wings of the Mgb lines, sensitive 
to changes in log jr. While we obtained good agreement with 
other log (/-sensitive spectral features (notably the wings of 
Ca I 6122 and Ca I 6162), as well as with the spectra of 
the Sun and Procyon, the synthetic profiles for the Mgb 
lines remained too wide. Determining \ogg from the Mg b 
lines alone yielded a very low value of logg = 3.95, which is 
strongly inconsistent with the photometrically determined 
radius of WASP-1. However, this apparent inconsistency 
can be reconciled by adjusting the elemental abundance of 
[Mg/H]. 

Because our selected wavelength regions do not contain 
Mg lines other than Mg b, we determined the Mg abun- 
dance in WASP-1 from seven weak Mg lines in other re- 
gions of the spectrum (Mg I 5528.4, 5711.1, 5785.3, 5785.6, 
6318.7, 6319.2 and 6319.5 A). We find [Mg/H]= 0.10, indi- 
cating that Mg is indeed underabundant by 0.12 in WASP-1. 
Similar underabundances in Mg (and Na) for metal rich 
stars and exoplanets hosts was reported by Gonzalez et al. 
(2001), although other comprehensive studies by Bensby et 
al. ( 2005 1 and Gilli et al. ( 2006 1 do not confirm such a trend. 



We then redetermined the stellar parameters of WASP-1, 
treating [Mg/H] as a fixed parameter. This produced only 
minor changes in the derived stellar parameters, but also a 
good agreement between the observed and synthetic profiles 
around the Mg b line complex. 

The final set of parameters we obtained from the analy- 
sis with SME are listed in Table [3~T| We also reproduce the 
synthetic spectra of the best choice of parameters in Fig [2] 
We would like to point out that all parameters were deter- 
mined simultaneously across all three wavelength regions, 
thus obtaining a self-consistent solution. 



3.3 Error estimates and parameter correlations 

As VF05 point out, the formal uncertainties reported by 
SME strongly underestimate the true uncertainties because 
the value of the goodness-of-fit parameter \ 2 is n °t domi- 
nated by the signal-to-noise of our spectra, but rather by 
uncertainties in the data reduction (for example, continuum 
fitting) as well as by inaccuracies in the synthetic spectra 
calculations. Thus, in analogy with VF05, we individually 
analysed our seven spectra of WASP-1 with SME, and quote 
the square root of the variance of the measurements of each 
parameter as the uncertainty in Table |3~T) 

Spectroscopically determined fundamental parameters 
tend to show strong correlations, and any uncertainty in 
one parameter may be propagated to other parameters. We 
investigated the magnitude of these correlations by varying 
T e ff with ±50 K, and determining a new optimal solution 
at these temperatures, thus propagating the uncertainty in 
temperature to the other parameters. The correlated uncer- 



tainties (listed in Table 3.1 1 indicate that the uncertainties 



out that it is unclear whether the assumption of LTE is applicable 
to the Balmer lines, and that non-LTE calculations might yield 
temperatures up to 100 K higher. 



due to parameter correlations are of similar order of magni- 
tude as the uncertainties derived earlier from the scatter of 
measurements of the seven individual spectra. 

In addition to propagating the uncertainties related to 
temperature, we also looked at the relation between [Mg/H], 
[M/H] and logp, as these may be expected to correlate 
through the Mg b line complex. We find a clear correla- 
tion between [Mg/H] and [M/H], such that [Mg/H] always 
shows a relative underabundance of about —0.1 with respect 
to [M/H]. On the other hand we see no indication that the 
value of log g correlates with [Mg/H] , mainly because each 
parameter is well constrained by diagnostics other than Mg 
b (i.e. the weak Mg lines and the broadening of Ca I 6122 
and Ca I 6162). 

One pair of parameters that shows a strong correlation 
are the assumed value of the macroturbulence v mac and the 
obtained projected rotational velocity i;sini. This is no sur- 
prise, because the observed photospheric line profiles are a 
convolution of rotational and turbulent broadening. In our 
analysis we have assumed a value of 4.5 kms -1 , similar to 
the choice of VF05. This value is in agreement with sev- 
eral other empirical relations for main-sequence type stars 
(see VF05). However, the parameters we derive for WASP- 
1 could indicate that WASP-1 is starting to evolve off the 
main-sequence (see Sect.[4]|, which would imply a slight in- 
crease in Wmac, to approximately 5.5 kms -1 . Using higher 
values for t> ma c would result in a lower value of the vsini, 
thus the value we quote in Table |3.1| should be considered 
an upper limit. 



3.4 Independent estimate of T e g from the IRFM 

In addition to the spectroscopic estimate, we also deter- 
mined the effective temperature using the IRFM. This 
method allows for a nearly model-independent determina- 
tion of T e g obtained from the integrated stellar flux at the 
Earth and a measurement at an infrared wavelength. 

For WASP-1 we used the Tycho B and V magnitudes 
and the 2MASS magnitudes to estimate the integrated flux 
(F 1 ®) by fitting Kurucz flux distributions. Since WASP-1 is 
relatively nearby and out of Galactic plane we assumed zero 
reddening. Integrating the fitted distributions gives a value 
of Ftot = (6.2 ± 0.5) ■ 10~ 13 W m" 2 . Using the IRFM to 
determine the T e g for the J, H and K s bands, we obtain 
an average T c g = 6200 ± 200 K, which is in close agreement 
with the spectroscopic determination. 



4 EVOLUTIONARY STATUS 

The stellar parameters derived from the detailed spec- 
troscopic analysis described above can be used to place 
WASP-1 in its evolutionary context, providing further con- 
straints on the age of the host star and its planet. Both the 
effective temperature and the (photometrically determined) 
stellar radius are measured to high accuracy using well un- 
derstood, robust physics that is largely model independent. 
Thus in Fig. [3] we plot a modified Hertzprung-Russel di- 
agram in which the radius acts as a proxy for luminosity 
(as it is related to luminosity through the temperature) and 
compare the properties of WASP-1 to several different super- 



solar metallicity stellar evolution models ( Girardi et al. 2002 
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Figure 3. Radius versus T e g (a modified Hertzprung-Russel di- 
agram) of the host stars of the 14 known transiting extra-solar 
planets [Burrows et al,| (asterisks, |2006^ , including WASP-1 (solid 
circle) . Here we present isochrones for super-solar metallicity stel- 
lar evolution model s from|Girardi et al.| ( |2002} (solid) , |Yi et ah] 
( |2001| > (dotted), and |Siess et al.| l |1997[ l (dashed). Given the uncer- 
tainties inherent in the models from various authors, the position 
of WASP-1 is consistent with an age of 1-3 Gyr. 



Siess et al.||1997| |Yi et al.||2001[ ) between 1-3 Gyr. The po- 
sition of WASP-1 in this diagram shows that its its age is 
likely to be in the range 1-3 Gyr and that its large stellar ra- 
dius may be explained by the fact that the star is beginning 
to evolve off the main sequence. 

We note that the mass initially adopted for WASP-1, 
M = 1.15 M Q , ( |Collier Cameron et aT1|2007[ ) is marginally 
inconsistent with the log g found here, when combined with 
the latest estimates of the stellar radius (i?* = 1.453 ± 



0.032i?. Charbonneau et al. 2007 Shporer et al. 20071. A 



slightly higher mass of M = 1.25-1.35 Mq would allow all 
the measured values of logg, temperature, and radius to be 
consistent given their errors. 

4.1 Lithium abundance 

An interesting feature in the spectrum of WASP-1 is the the 
strong lithium absorption line at 6708 A. As stars only de- 
stroy lithium, the amount of photospheric lithium is often 
used as an age indicator. However, the accuracy of age es- 
timates based on lithium abundances depends strongly on 
the efficiency with which stars deplete lithium. Stars such 
as WASP-1, with temperatures of ~ 6100 K, cannot deplete 
lithium very easily from their surface layers, because the 
bottom of the convection zone does not reach temperatures 
high enough to burn lithium. It is also cooler than 6300 K, 
the lower limit of the so-called "Li dip" or "Boesgaard gap" 



(for example, Balachandran 1995 1 where lithium depletion 
is enhanced by an as yet unknown process. Still, early type 
stars do slowly deplete lithium during their main sequence 
lifetimes ( Jones et al.|1999 l, allowing modest constraints to 
be placed on the age on WASP-1 based on the abundance 
of lithium present in the photosphere. 

In the spectrum of WASP-1, Li I 6708 has an equiva- 
lenth width of 135 mA. Using SME and the atmospheric 
parameters derived earlier, we find that this corresponds 




log n(Li) = 2.91 ± 0.05 



6707.0 6707.5 6708.0 6708.5 

Wavelength (A) 

Figure 4. A range of synthetic line profiles for the Li I 6708 
feature, in steps of 0.05 in logn(Li). The thick line corresponds 
to the optimum value, logn(Li) = 2.91. Changes of ±50 K in 
the stellar effective temperature result in changes of ±0.05 in the 
lithium abundance (see text). 



to an elemental abundance of (Li/H) = —9.09 ± 0.05, or 
logn(Li) = (Li/H) + 12 = 2.91 ± 0.05 in more conventional 
units (see Figure Eh. We find a weak correlation between 
T e ff and logn(Li), where a decrease in temperature by 50 K 
relates to a decrease in logn(Li) by 0.05. 

The value logn(Li) = 2.91 ± 0.05 is less than the com- 
monly adopted value of 3.2 for the primordial lithium abun- 
dance. Still, it is a the high end with respect to the study by 
Israelian et al. 1(12004 



who compared the amount of lithium 
in exoplanet host stars with a set of reference stars. We also 
compared our value of logn(Li) with the empirical relations 
of lithium depiction in open clusters derived by Sestito fc 
Randich ( |2005 I. Their study shows that, for stars with tem- 
peratures of 6200 ± 150 K, lithium can be expected to reach 
logn(Li) = 2.9 after 1-3 Gyr. This age estimate is consistent 
with what we derived earlier from evolutionary tracks. 



5 DISCUSSION 

Our spectroscopic analysis, together with the stellar mass 
and radius estimates inferred from previous studies, reveal 
a remarkable similarity between the physical parameters of 
WASP-1 and HD 149026 (see Table [5] The two stars appear 
to be very similar in both mass and evolutionary status. 



Sato et al. (20051 estimate the age of HD 149026 to be 2.0 ± 



0.8 Gyr. Also the environments in which the two planets 
orbit are remarkably similar. WASP- lb orbits slightly closer 
to a slightly cooler star, with the result that its blackbody 
effective temperature (assuming a low albedo and efficient 
redistribution and reradiation of heat from the night side) 
is just 80 K warmer than that of HD 149026b. 

The structural differences between the two planets con- 
trasts starkly with the similarities of the parent stars. The 
light, Saturn-mass HD 149026b requires a core mass of sev- 
eral tens of Earth masses to reproduce its small radius and 
high density ( |Sato et al 1 2005)1 Fort ney et al.|2006| ). WASP- 
lb is more than twice as massive and has a much larger 
radius. This can only be reconciled with the models of |Bur-| 
rows et al. ( 2006 1 if it has little or no core, its atmosphere 



is metal-rich, and its age is less than 1.5 Gyr. 

Although HD 149026 could be a factor 2 more metal- 
rich than WASP-1 within the uncertainties, WASP-1 itself 
is substantially more metal-rich than the Sun. With this 
caveat in mind, we suggest that these two planets provide 
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Table 2. A comparison between the stellar and planetary param- 
eters ofWASP-l_andHD 14902 6. D ata for HP 149026 are taken 
from |Charbonneau et ai7| (|2006[l and|Sato et al.| ( |2005[l. Data for 
WASP-1 are from |ColIier Cameron et al.| ( |2007[ i"and |Charbonneau| 
|et al.| ( [2007| , except T e ff )t , [Fe/H] and the age estimate, which are 
from this paper. The values of T e ff p were determined assuming 
zero albedo and full redistribution of heat to the nightside. 



HD149026 



WASP-1 





(M@) 


1.3 ± 0.1 


i 

L - 10 -0.09 


R, 


(Re) 


1.45 ± 0.1 


1.453 ± 0.032 




(K) 


6147 ± 50 K 


6110 ± 45 K 


[Fe/H] 




0.36 ± 0.05 


0.26 ± 0.03 


Age 


(Gyr) 


2.0 ± 0.8 


2.0 ± 1.0 


M p 


(Mj up ) 


0.330 ± 0.002 


7Q+013 

u - ' a -0.06 


Rp 


(^?Jup) 


0.726 ± 0.064 


1.443 ± 0.039 


Pv 


(PJup) 


0.86 


0.26 


orbital sep. 


(AU) 


0.042 


0.0379 ± 0.0042 


T e ft,p 


(K) 


1740 


1820 



a strong counter-example to the trend suggested tentatively 
by |GuiUot et al.| ( |2006] ) and |Burrows et al.| ( [2006] ), in which 
core mass increases with metallicity of the host star. We 



conclude that the final balance between core and envelope 
mass in giant planets cannot be a simple function of stellar 
metallicity, nor of radiation environment. As |Ikoma et al.| 
(2006) conclude, the final composition of a planet is more 
likely to be dictated by the details of the disk environment in 
which the planet formed, and possibly the dynamical history 
following the envelope accretion phase, than by the compo- 
sition of the parent star. 



ACKNOWLEDGMENTS 

We would like to thank the referees for constructive criti- 
sism that helped to improve this paper. HCS acknowledges 
support from the Swedish Research Council. The Nordic Op- 
tical Telescope is operated on the island of La Palma jointly 
by Denmark, Finland, Iceland, Norway, and Sweden, in the 
Spanish Observatorio del Roque de los Muchachos of the 
Instituto de Astrofisica de Canarias. 



REFERENCES 

Bakos G. A., et al., 2007, ApJ, 656, 552 
Balachandran S., 1995, ApJ, 446, 203 

Barklem P. S., 2007, A&A, in press, arXiv:astro- 
|ph/0702222] 

Barklem P. S., Piskunov N., O'Mara B. J., 2000, A&A, 363, 
1091 

Barklem P. S., Stempels H. O, Allende Prieto O, 
Kochukhov O. P., Piskunov N., O'Mara B. J., 2002, A&A, 
385, 951 

Bensby T., Feltzing S., Lundstrom I., Ilyin L, 2005, A&A, 
433, 185 

Blackwell D. E., Petford A. D., Shallis M. J., 1980, A&A, 
82, 249 

Bodenheimer P., Lin D. N. O, Mardling R. A., 2001, ApJ, 
548, 466 



Burrows A., Hubeny L, Bu daj J., Hubbar d W. B., 2006, 

ApJ, submitted, |arXiv:astro- ph/0612703"| 
Charbonneau D., Brown T. M., Latham D. W., Mayor M., 

2000, ApJ, 529, L45 
Charbonneau D., et al., 2006, ApJ, 636, 445 
Charbonneau D., Winn J. N., Everett M. E., Latham 

D. W., Holman M. J., Esquerdo G. A., O'Donovan F. T., 

2007, ApJ, 658, 1322 
Collier Cameron A., et al., 2007, MNRAS, 375, 951 
Fortney J. J., Saumon D., Marley M. S., Lodders K., Freed- 

man R. S., 2006, ApJ, 642, 495 
Fuhrmann K., Pfeiffer M., Frank O, Reetz J., Gehren T., 

1997, A&A, 323, 909 
Girardi L., Bertelli C, Bressan A., Chiosi O, Groenewegen 

M. A. T., Marigo P., Salasnich B., Weiss A., 2002, A&A, 

391, 195 

Gilli C, Israelian G., Ecuvillon A., Santos N. O, Mayor 

M., 2006, A&A, 449, 723 
Gonzalez G., Laws O, Tyagi S., Reddy B. E., 2001, AJ, 

121, 432 

Guillot T., Showman A. P., 2002, A&A, 385, 156 

Guillot T., Santos N. C, Pont F., Iro N., Melo O, Ribas 

L, 2006, A&A, 453, L21 
Henry G. W., Marcy G. W., Butler R. P., Vogt S. S., 2000, 

ApJ, 529, L41 

Ikoma M., Guillot T., Genda H., Tanigawa T., Ida S., 2006, 

ApJ, 650, 1150 
Israelian G., Santos N. O, Mayor M., Rebolo R., 2004, 

A&A, 414, 601 
Jones, B. F., Fischer, D., & Soderblom, D. R. 1999, AJ, 

117, 330 

Kupka F., Piskunov N., Ryabchikova T. A., Stempels H. O, 

Weiss W. W., 1999, A&AS, 138, 119 
Kurucz R., 1993, ATLAS9 Stellar Atmosphere Programs 

and 2 km/s grid. (Kurucz CD-ROM No. 13) Cambridge, 

Mass.: Smithsonian Astrophysical Observatory 
Kurucz R. L., Furenlid I., Brault J., Testerman L., 1984, 

NSO Atlas, Sunspot, New Mexico: National Solar Obs. 
McCullough P. R., et al., 2006, ApJ, 648, 1228 
Piskunov N. E., Kupka F., Ryabchikova T. A., Weiss 

W. W., Jeffery C. S., 1995, A&AS, 112, 525 
Santos N. C, Israelian C, Mayor M., 2004, A&A, 415, 1153 
Santos N. O, et al., 2006, A&A, 458, 997 
Sato B., et al., 2005, ApJ, 633, 465 
Sestito P., Randich S., 2005, A&A, 442, 615 
Shporer A., Tamuz P., Zucker S., M azeh T., 2007, MNRAS, 

in press, arXiv:astro-ph/0610556 
Siess, L., Forestini, M., & Dougados, C. 1997, A&A, 324, 

556 

Valenti J. A., Fischer D. A., 2005, ApJS, 159, 141 
Valenti J. A., Piskunov N., 1996, A&AS, 118, 595 
Winn J. N., Holman M. J., 2005, ApJ, 628, L159 
Yi, S., Demarque, P., Kim, Y., Lee, Y., Ree, C. H., Lejeune, 
T., & Barnes, S. 2001, ApJS, 136, 417 

This paper has been typeset from a TjnX/ DTjtX file prepared 
by the author. 



